Diffractive Dijet Production at HERA 
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Abstract 

We present recent experimental data from the HI and ZEUS Collaborations 
at HERA for diffractive dijet production in deep-inelastic scattering (DIS) and 
photoproduction and compare them with next-to-leading order (NLO) QCD 
predictions using diffractive partem densities. While good agreement is found 
for DIS, the dijet photoproduction data are overestimated by the NLO theory, 
showing that factorization breaking occurs at this order. While this is expected 
theoretically for resolved photoproduction, the fact that the data are better de- 
scribed by a global suppression of direct and resolved contribution by about 
a factor of two comes as a surprise. We therefore discuss in some detail the 
factorization scheme and scale dependence between direct and resolved con- 
tributions and propose a new factorization scheme for diffractive dijet photo- 
production. 

1 Introduction 

It is well known that in high-energy deep-inelastic ep-collisions a large fraction of the observed events are 
diffractive. These events are defined experimentally by the presence of a forward-going system Y with 
four-momentum py, low mass My (in most cases a single proton and/or low-lying nucleon resonances), 
small momentum transfer squared t = {p — Py) 2 , and small longitudinal momentum transfer fraction 
X P = lip ~ Py)/qp from the incoming proton with four-momentum p to the system X (see Fig.Q. 
The presence of a hard scale, as for example the photon virtuality Q 2 = —q 2 in deep-inelastic scattering 
(DIS) or the large transverse jet momentum p* T in the photon-proton centre-of-momentum frame, should 
then allow for calculations of the production cross section for the central system X with the known meth- 
ods of perturbative QCD. Under this assumption, the cross section for the inclusive production of two 
jets, e + p — > e + 2 jets + X' + Y, can be predicted from the well-known formulae for jet production in 




Fig. 1: Diffractive scattering process ep — > eXY, where the hadronic systems X and Y are separated by the largest rapidity 
gap in the final state. 



non-diffractive ep collisions, where in the convolution of the partonic cross section with the partem distri- 
bution functions (PDFs) of the proton the latter ones are replaced by the diffractive PDFs. In the simplest 
approximation, they are described by the exchange of a single, factorizable pomeron/Regge-pole. 

The diffractive PDFs have been determined by the HI Collaboration at HERA from high-precision 
inclusive measurements of the DIS process ep — > eXY using the usual DGLAP evolution equations in 
leading order (LO) and next-to-leading order (NLO) and the well-known formula for the inclusive cross 
section as a convolution of the inclusive parton-level cross section with the diffractive PDFs [1]. For a 
similar analysis of the inclusive measurements of the ZEUS Collaboration see El|3]]. A longer discussion 
of the extraction of diffractive PDFs can also be found in these proceedings [4] and in [5]. For inclusive 
diffractive DIS it has been proven by Collins that the formula referred to above is applicable without 
additional corrections and that the inclusive jet production cross section for large Q 2 can be calculated 
in terms of the same diffractive PDFs [6]. The proof of this factorization formula, usually referred to as 
the validity of QCD factorization in hard diffraction, may be expected to hold for the direct part of pho- 
toproduction (Q 2 ~ 0) or low-Q 2 electroproduction of jets |6|. However, factorization does not hold for 
hard processes in diffractive hadron-hadron scattering. The problem is that soft interactions between the 
ingoing two hadrons and their remnants occur in both the initial and final state. This agrees with experi- 
mental measurements at the Tevatron [7 ]. Predictions of diffractive dijet cross sections for pp collisions 
as measured by CDF using the same PDFs as determined by HI [1| overestimate the measured cross 
section by up to an order of magnitude [7 1. This suppression of the CDF cross section can be explained 
by considering the rescattering of the two incoming hadron beams which, by creating additional hadrons, 
destroy the rapidity gap [ 8 ]. 

Processes with real photons (Q 2 ~ 0) or virtual photons with fixed, but low Q 2 involve direct 
interactions of the photon with quarks from the proton as well as resolved photon contributions, leading 
to parton-parton interactions and an additional remnant jet coming from the photon (for a review see 
1191). As already said, factorization should be valid for direct interactions as in the case of DIS, whereas 
it is expected to fail for the resolved process similar as in the hadron-hadron scattering process. In a 
two-channel eikonal model similar to the one used to calculate the suppression factor in hadron-hadron 
processes [8], introducing vector-meson dominated photon fluctuations, a suppression by about a factor 
of three for resolved photoproduction at HERA is predicted [ 10 1. Such a suppression factor has recently 
been applied to diffractive dijet photoproduction fTTllT2l and compared to preliminary data from HI f\3\ 
and ZEUS [ 14 1. While at LO no suppression of the resolved contribution seemed to be necessary, the 
NLO corrections increase the cross section significantly, showing that factorization breaking occurs at 
this order at least for resolved photoproduction and that a suppression factor R must be applied to give a 
reasonable description of the experimental data. 

As already mentioned elsewhere fTTllL? !. describing the factorization breaking in hard photopro- 
duction as well as in electroproduction at very low Q 2 lfT31 by suppressing the resolved contribution 
only may be problematic. An indication for this is the fact that the separation between the direct and the 
resolved process is uniquely defined only in LO. In NLO these two processes are related. The separation 
depends on the factorization scheme and the factorization scale M 7 . The sum of both cross sections 
is the only physically relevant cross section, which is approximately independent of the factorization 
scheme and scale [ 16 1. As demonstrated in Refs. iTTTl IT2l multiplying the resolved cross section with 
the suppression factor R = 0.34 destroys the correlation of the M 7 -dependence between the direct and 
resolved part, and the sum of both parts has a stronger M 7 -dependence than for the unsuppressed case 
(R = 1), where the M 7 -dependence of the NLO direct cross section is compensated to a high degree 
against the M 7 -dependence of the LO resolved part. 



In the second Section of this contribution, we present the current experimental data from the HI 
and ZEUS Collaborations on diffractive dijet production in DIS and photoproduction and compare these 
data to theoretical predictions at NLO for two different scenarios: suppression of only the resolved part 
by a factor R = 0.34 as expected from LO theory and proposed in [ 8 1, and equal suppression of all direct 
and resolved contributions by a factor R = 0.5, which appears to describe the data better phenomeno- 
logically. This motivates us to investigate in the third Section the question whether certain parts of the 
direct contribution might break factorization as well and therefore need a suppression factor. 

The introduction of the resolved cross section is dictated by perturbation theory. At NLO, collinear 
singularities arise from the photon initial state, which are absorbed at the factorization scale into the pho- 
ton PDFs. This way the photon PDFs become Af 7 -dependent. The equivalent M 7 -dependence, just with 
the opposite sign, is left in the NLO corrections to the direct contribution. With this knowledge, it is 
obvious that we can obtain a physical cross section at NLO, i.e. the superposition of the NLO direct 
and LO resolved cross section, with a suppression factor R < 1 and no M 7 -dependence left, if we also 
multiply the In M 7 -dependent term of the NLO correction to the direct contribution with the same sup- 
pression factor as the resolved cross section. We are thus led to the theoretical conclusion that, contrary 
to what one may expect, not all parts of the direct contribution factorize. Instead, the initial state singular 
part appearing beyond LO breaks factorization even in direct photoproduction, presumably through soft 
gluon attachments between the proton and the collinear quark-antiquark pair emerging from the photon 
splitting. This would be in agreement with the general remarks about initial state singularities in Ref. 0. 

In the third Section of this contribution, we present the special form of the In M 7 -term in the NLO 
direct contribution and demonstrate that the M 7 -dependence of the physical cross section cancels to a 
large extent in the same way as in the unsuppressed case (R = 1). These studies can be done for pho- 
toproduction (Q 2 ~ 0) as well as for electroproduction with fixed, small Q 2 . Since in electroproduction 
the initial-state singularity in the limit Q 2 — > is more directly apparent than for the photoproduction 
case, we shall consider in this contribution the low-Q 2 electroproduction case just for demonstration. 
This diffractive dijet cross section has been calculated recently JT31 . A consistent factorization scheme 
for low-Q 2 virtual photoproduction has been defined and the full (direct and resolved) NLO corrections 
for inclusive dijet production have been calculated in 1171 . In this work we adapt this inclusive NLO 
calculational framework to diffractive dijet production at low-Q 2 in the same way as in fi"5ll . except that 
we multiply the In M 7 -dependent terms as well as the resolved contributions with the same suppression 
factor R = 0.34, as an example, as in our earlier work fTTl [1211131. The exact value of this suppression 
factor may change in the future, when better data for photoproduction and low-Q 2 electroproduction 
have been analyzed. We present the In M 7 -dependence of the partly suppressed NLO direct and the fully 
suppressed NLO resolved cross section da / dQ 2 and their sum for the lowest Q 2 bin, before we give a 
short summary in section 4. 

2 Comparison of HI and ZEUS Data with NLO Theory Predictions 

In this Section, diffractive PDFs Q |2j Q extracted from diffractive structure function data are used in 
NLO calculations to test factorisation in diffractive dijet production. Dijet production is directly sensitive 
to the diffractive gluon (Fig. |2]) whereas in inclusive measurements the gluon is determined from scaling 
violations. 



2.1 Diffractive Dijet Production in DIS 

HI has measured the cross sections for dijet pro- 
duction [ 13 1 in the kinematic range Q 2 > 4 GeV 2 , 
165 < W < 242 GeV (photon-proton centre-of- 
mass energy) and xjp < 0.03. Jets are identified 
using the inclusive kx cluster algorithm and se- 
lected by requiring £$ jet (l, 2) > 5,4 GeV and 
—3 < r/j* ct < 0. 1 NLO predictions have been ob- 
tained by interfacing the HI diffractive PDFs with 
the DISENT program [18]. The renormalisation 
and factorisation scales were set to the transverse 
energy of the leading parton jet. The NLO parton jet cross sections have been corrected for hadronisation 
effects using the leading order (LO) generator RAPGAP [ 19] with parton showers and the Lund string 
fragmentation model. Comparisons of the DISENT and RAPGAP predictions with the measured cross 
section differential in 2p ts , an estimator for the fraction of the momentum of the diffractive exchange 
entering the hard scatter, are shown in Fig. |3^. The inner band around the NLO calculation indicates 





Fig. 2: Example processes for a) direct photon and b) resolved 
photon interactions. 
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Fig. 3: Diffractive DIS dijet cross sections compared with a NLO prediction based on diffractive PDFs and with RAPGAP. 



the « 20% uncertainty resulting from a variation of the renormalisation scale by factors 0.5 and 2. The 
uncertainty in the diffractive PDFs is not shown. Within this additional uncertainty, which is large at 
high z^ s , the cross section is well described. The cross section differential in log 10 (xjp), p 1 ^ 1 , and Q 2 
is shown in Fig. fJJ) and Fig. |4] All distributions are well described and QCD factorisation is therefore in 
good agreement with dijet production in diffractive DIS. 

Similar results are presented by ZEUS [20|; the dijet cross sections have been measured in the 
kinematic range 5 < Q 2 < 100 GeV 2 , 100 < W < 200 GeV, xjp < 0.03. The jets were identified 
using the inclusive kx algorithm in the 7p frame and required to satisfy E T ' ie (1, 2) > 5,4 GeV and 
— 3.5 < 77j* ct < 0.0. NLO predictions have been obtained with the DISENT program interfaced to three 
different sets of diffractive PDFs: from fit to HI data Q, from fit to the ZEUS M x data (GLP) 
and from fit to ZEUS LPS and p^ charm d a t a jjj. Comparisons of the DISENT predictions with the 
measured cross section differential in E T ' ]et , r]? et , Zp s and x° bs are shown in Fig. |5] The 20 — 30% 
uncertainty in the NLO calculations resulting from a variation of the renormalisation and factorisation 
scales is not shown. Within the experimental and QCD scale uncertainties, the predictions based on 
the HI and ZEUS-LPS PDFs give a good description of the dijet cross section. The normalisation of 

'The '*' denotes variables in the photon-proton centre-of-mass system. 
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Fig. 4: Diffractive DIS dijet cross sections compared with a NLO prediction based on diffractive PDFs. 



the prediction using the GLP tit is substantially lower than those from the other two sets of PDFs. For 
ZEUS, the difference observed between the three sets may be interpreted as an estimate of the uncertainty 
associated with the diffractive PDFs and with the definition of the diffractive region. The dijet data could 
be included in future fits in order to better constrain the diffractive gluon distribution. 

Within the experimental and theoretical uncertainties and assuming the HI diffractive PDFs, fac- 
torisation is in good agreement with diffractive D* production [21 22 1 in the DIS kinematic region. 

2.2 Diffractive Photoproduction of Dijets 

In photoproduction, a sizeable contribution to the cross section is given by resolved photon processes 
(Fig. I2J5) in which only a fraction x 7 < 1 of the photon momentum enters the hard scatter. The photo- 
production dijet cross section measured by HI (Q 2 < 0.01 GeV 2 , 165 < W < 242 GeV, x P < 0.03, 
(1, 2) > 5, 4 GeV, — 1 < r/j ct < 2, inclusive kx algorithm) is shown in Fig.|6][13]. NLO predictions 
have been obtained with the Frixione et al. program! 23 1 interfaced to the HI diffractive PDFs. The 
parton jet calculation is corrected for hadronisation effects using RAPGAP. The cross section differential 
in 2Tp ts is shown in Fig. [6^. The calculation lies a factor rs 2 above the data. Fig. |6p and|6j; show the 
cross section as a function of p^ 1 and x J 7 ets and the NLO predictions have been scaled down by a factor 
0.5. Good agreement is obtained for this global suppression. In Fig. |6jl, only the "resolved" part for 
which x J 7 ets < 0.9 at the parton level is scaled by a factor 0.34. This factor was proposed by Kaidalov et 
al. CO I for the suppression of the resolved part in LO calculations. The calculation for x^ ts > 0.9 is left 
unsealed. This approach is clearly disfavoured. 

The ZEUS measurement |24) (Q 2 < 0.01 GeV 2 , X]P < 0.025, 0.2 < y < 0.85, £^(1,2) > 
7.5,6.5 GeV, —1.5 < 77 < 1.5, inclusive kx algorithm) is shown in Fig. and Fig. |8] separately for 
samples enriched in "direct" (x J 7 ets > 0.75) and "resolved" (x J 7 ets < 0.75) processes, respectively. The 
NLO [ 12 1 prediction using the HI diffractive PDFs is also presented corrected for hadronization effects 
and with the "resolved" part scaled by the factor 0.34. No evidence is observed for a suppression of 
resolved photon processes relative to direct photon processes in any particular kinematic region. 

Diffractive dijet photoproduction is overestimated by calculations based on PDFs which give a 
good description of the diffractive DIS data. Factorisation is broken in photoproduction relative to DIS 
by a factor m 0.5 with no observed dependence on x 7 or other kinematic variables. 
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Fig. 5: Diffractive DIS dijet cross sections compared with NLO predictions based on three sets of diffractive PDFs. 



3 Factorization and its Breaking in Diffractive Dijet Production 

The fact that equal suppression of direct and resolved photoproduction by a factor R = 0.5 appears to 
describe the HI and ZEUS data better phenomenologically motivates us to investigate in some detail the 
question whether certain parts of the direct contribution might break factorization as well and therefore 
need a suppression factor. These studies can be done for photoproduction (Q 2 ~ 0) as well as for elec- 
troproduction with fixed, small Q 2 . Since in electroproduction the initial-state singularity in the limit 
Q 2 —* is more directly apparent than for the photoproduction case, we shall consider in this contribu- 
tion the low-Q 2 electroproduction case just for demonstration. 



A factorization scheme for virtual photoproduction has been defined and the full NLO corrections 
for inclusive dijet production have been calculated in lUTl. They have been implemented in the NLO 
Monte Carlo program JETVIP [25 1 and adapted to diffractive dijet production in fT"5l. The subtraction 
term, which is absorbed into the PDFs of the virtual photon / a / 7 (x 7 , M 7 ), can be found in \26\. The 
main term is proportional to ln(M 2 / Q 2 ) times the splitting function 

P q ^(z) = 2N c Q 2 z2 + { \- z) \ (1) 

where z = P1P2/P0Q £ [%', 1] and Qi is the fractional charge of the quark q^. p\ and P2 are the momenta 
of the two outgoing jets, and po and q are the momenta of the ingoing partem and virtual photon, respec- 
tively. Since Q 2 = —q 2 <C M 2 , the subtraction term is large and is therefore resummed by the DGLAP 
evolution equations for the virtual photon PDFs. After this subtraction, the finite term M(Q 2 )y^, which 
remains in the matrix element for the NLO correction to the direct process fTTl. has the same in- 
dependence as the subtraction term, i.e. lnM 7 is multiplied with the same factor. As already mentioned, 
this yields the M 7 -dependence before the evolution is turned on. In the usual non-diffractive dijet photo- 
production these two M 7 -dependences cancel, when the NLO correction to the direct part is added to the 
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Fig. 6: a) Diffractive dijet photoproduction cross section differential in z^, ts compared with a NLO prediction based on diffrac- 
tive PDFs and RAPGAP. b)-d): Cross section differential in p J y U and x^ ts , compared with the NLO prediction modified as 
follows: in b) and c) the calculation is scaled by a global factor 0.5 whereas in d) only the "resolved" part is scaled by 0.34. 



LO resolved cross section HI 61 . Then it is obvious that the approximate Af 7 -independence is destroyed, if 
the resolved cross section is multiplied by a suppression factor R to account for the factorization break- 
ing in the experimental data. To remedy this deficiency, we propose to multiply the In M 7 -dependent 
term in M(Q 2 )y^ with the same suppression factor as the resolved cross section. This is done in the 
following way: we split M(Q 2 )j^g into two terms using the scale p^ in such a way that the term contain- 
ing the slicing parameter y s , which was used to separate the initial-state singular contribution, remains 
unsuppressed. In particular, we replace the finite term after the subtraction by 
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where R is the suppression factor. This expression coincides with the finite term after subtraction (see 
Ref. [ 26 1) for R = 1, as it should, and leaves the second term in Eq. © unsuppressed. In Eq. (0 we have 
suppressed in addition to ln(M 2 /p^) also the z-dependent term ln(z/(l — z)), which is specific to the 
MS subtraction scheme as defined in fTD. The second term in Eq. © must be left in its original form, 
i.e. being unsuppressed, in order to achieve the cancellation of the slicing parameter (y s ) dependence 
of the complete NLO correction in the limit of very small Q 2 or equivalently very large s. It is clear 
that the suppression of this part of the NLO correction to the direct cross section will change the full 
cross section only very little as long as we choose M 7 ~ pfp. The first term in Eq. which has the 
suppression factor R, will be denoted by DIRis in the following. 
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Fig. 7: Direct enriched photoproduction. Diffractive dijet photoproduction cross section differential in y, xp, z 3 p S , 

and 77j ct x compared with a NLO prediction based on diffractive PDFs. The NLO prediction is also presented corrected for 

hadronization effects and with the "resolved" part scaled by 0.34. 



To study the left-over M 7 -dependence of the physical cross section, we have calculated the diffrac- 
tive dijet cross section with the same kinematic constraints as in the HI experiment [27 ]. Jets are defined 
by the CDF cone algorithm with jet radius equal to one and asymmetric cuts for the transverse momenta 
of the two jets required for infrared stable comparisons with the NLO calculations 1281 . The original HI 
analysis actually used a symmetric cut of 4 GeV on the transverse momenta of both jets [29 1. The data 
have, however, been reanalyzed for asymmetric cuts [27 1. 

For the NLO resolved virtual photon predictions, we have used the PDFs SaSID [30] and trans- 
formed them from the DIS 7 to the MS scheme as in Ref. fPTI . If not stated otherwise, the renormalization 
and factorization scales at the pomeron and the photon vertex are equal and fixed to pfj, = p^ - etx . We 
include four flavors, i.e. nj = 4 in the formula for a s and in the PDFs of the pomeron and the photon. 
With these assumptions we have calculated the same cross section as in our previous work [ 15 1. First we 
investigated how the cross section da/dQ 2 depends on the factorization scheme of the PDFs for the vir- 
tual photon, i.e. da/dQ 2 is calculated for the choice SaSID and SaSlM. Here da/dQ 2 is the full cross 
section (sum of direct and resolved) integrated over the momentum and rapidity ranges as in the HI anal- 
ysis. The results, shown in Fig. 2 of Ref. [26], demonstrate that the choice of the factorization scheme 
of the virtual photon PDFs has negligible influence on da/dQ 2 for all considered Q 2 . The predictions 
agree reasonably well with the preliminary HI data [.2*71 . 

We now turn to the M 7 -dependence of the cross section with a suppression factor for DIRis. To 
show this dependence for the two suppression mechanisms, (i) suppression of the resolved cross section 
only and (ii) additional suppression of the DIRig term as defined in Eq. © in the NLO correction of 
the direct cross section, we consider da/dQ 2 for the lowest Q 2 -bin, Q 2 € [4,6] GeV 2 . In Fig. [9j 
this cross section is plotted as a function of £ = M^/p^ in the range £ 6 [0.25; 4] for the cases (i) 
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Fig. 8: Resolved enriched photoproduction. Diffractive dijet photoproduction cross section differential in y, xp, 2p S , 

and 77j ct 1 compared with a NLO prediction based on diffractive PDFs. The NLO prediction is also presented corrected for 

hadronization effects and with the "resolved" part scaled by 0.34. 



(light full curve) and (ii) (full curve). We see that the cross section for case (i) has an appreciable 
^-dependence in the considered £ range of the order of 40%, which is caused by the suppression of 
the resolved contribution only. With the additional suppression of the DIRis term in the direct NLO 
correction, the ^-dependence of da/dQ 2 is reduced to approximately less than 20%, if we compare the 
maximal and the minimal value of da/dQ 2 in the considered £ range. The remaining ^-dependence 
is caused by the NLO corrections to the suppressed resolved cross section and the evolution of the 
virtual photon PDFs. How the compensation of the M 7 -dependence between the suppressed resolved 
contribution and the suppressed direct NLO term works in detail is exhibited by the dotted and dashed- 
dotted curves in Fig. |9] The suppressed resolved term increases and the suppressed direct NLO term 
decreases by approximately the same amount with increasing £. In addition we show also da/dQ 2 in the 
DIS theory, i.e. without subtraction of any InQ 2 terms (dashed line). Of course, this cross section must 
be independent of £. This prediction agrees very well with the experimental point, whereas the result for 
the subtracted and suppressed theory (full curve) lies slightly below. We notice, that for M 7 = the 
additional suppression of DIRis has only a small effect. It increases da/dQ 2 by 5% only. 



4 Summary 

Experimental data from the HI and ZEUS Collaborations at HERA for diffractive dijet production in DIS 
and photoproduction have been compared with NLO QCD predictions using diffractive parton densities 
from HI and ZEUS. While good agreement was found for DIS assuming the HI diffractive PDFs, the 
dijet photoproduction data are overestimated by the NLO theory, showing that factorization breaking 
occurs at this order. While this is expected theoretically for resolved photoproduction, the fact that the 
data are better described by a global suppression of direct and resolved contribution by about a factor 
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Fig. 9: Photon factorization scale dependence of resolved and direct contributions to da/dQ 2 together with their weighted 
sums for (i) suppression of the resolved cross section and for (ii) additional suppression of DIRis, using SaSID virtual photon 
PDFs (30|. 

of two has come as a surprise. We have therefore discussed in some detail the factorization scheme and 
scale dependence between direct and resolved contributions and proposed a new factorization scheme 
for diffractive dijet photoproduction. 
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